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Abstract 

We study the two-dimensional periodic Anderson model at half-filling us- 
ing quantum Monte Carlo (QMC) techniques. The ground state undergoes a 
magnetic order-disorder transition as a function of the effective exchange cou- 
pling between the conduction and localized bands. Low-lying spin and charge 
excitations are determined using the maximum entropy method to analyti- 
cally continue the QMC data. At finite temperature we find a competition 
between the Kondo effect and antiferromagnetic order which develops in the 

localized band through Ruderman-Kittel-Kasuya-Yosida interactions. 
PACS numbers: 75.20.Hr, 75.30.Et, 75.30.Mb, 75.40.Mg 
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The periodic Anderson model (PAM) |T| describes a localized /-band of strongly cor- 
related electrons hybridized with a d-band of free conduction electrons. As the tempera- 
ture, filling, and Hamiltonian parameters are varied, there is a competition between the 
Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction [[| and the Kondo effect ||. The 
RKKY interaction favors ordering the magnetic moments of the localized /-band, while the 
Kondo effect screens the localized magnetic moments and quenches the magnetic interaction 
through the formation of singlets between the two bands [|J. The Kondo effect is expected 
to dominate for large exchange coupling between the conduction electrons and the local 
moments. For smaller values of the exchange interaction, instead, the local moments order, 
in most cases antiferromagnetically. 

It is believed that this behavior qualitatively describes the competition between magnetic 
ordering and singlet formation in a number of the heavy fermion materials. Additionally, 
several small-gap semiconducting compounds involving either rare-earth or actinide metals 
H can be well described by the PAM. These systems exhibit local moment behavior at 
high temperatures. At low temperatures the hybridization of the local moments with the 
conduction band leads to magnetic ordering and the formation of a small energy gap or, in 
some cases, to a Kondo insulator with no long-range magnetic order. 

The competition between the RKKY interaction and the Kondo effect has been studied 
within the framework of the two- impurity Anderson model P-flOl , the one-dimensional PAM 



TT| and the Kondo lattice ||12|| , and in infinite dimensions fl3|| . Here we present QMC results 
for the two-dimensional PAM at half-filling. This is an interesting model since it exhibits 
various types of insulating states ranging from an antiferromagnetic insulator to a Kondo 
insulator or, if the f-d hybridization dominates, a simple band- insulator. 
The Hamiltonian for the two-dimensional PAM is 



H = -tJ2 {4 a d ja + d)j ia ) + e d £ nt + U f £ 

(i,j)a iv i 



/ l \ Lf 1 



+e f £ n{ CT - V £ {4j ia + fUia) ■ (1) 
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Here t is the hopping parameter in the d-band, V is the hybridization energy between the 



bands, Uf is the Coulomb repulsion in the /-band, e/ and td are the energy levels of the el- 
and /-band respectively, and nf a = d\ a d ia and n{ a = fj a fi a are the density operators for the 
two bands. In the following we will set t — 1 and consider the symmetric case (e/ = 0) at 
half-filling (e d = 0). 

The zero temperature behavior of the PAM along the Uf — and V = axes in the U f-V 
phase diagram is simple. In the former case, we have a "band-insulator" as the hybridization 
V opens up a gap where the dispersionless / electron level crosses the conduction band. In 
the latter case we have a set of individual local moments completely decoupled from the 
conduction band. When V and U / are both non-zero, a qualitative understanding of the 
zero temperature order- disorder transition can be obtained from a comparison of the RKKY 
energy scale J RKK y ~ J 2 /W and the Kondo energy scale T K ~ We~ w / J , where J ~ V 2 /Uf 
and W is the bandwidth. For small values of J, long-range RKKY driven antiferromagnetic 
order will set in. As the exchange coupling J is increased beyond a fixed value, we expect 
that the order will be quickly suppressed due to the formation of spin singlets between the 
two bands. The above argument suggests a functional form of the phase boundary with 
U f oc V 2 . 

In order to determine whether the ground state of the PAM does have long-range an- 
tiferromagnetic order, we have carried out a finite size scaling analysis of the equal time 
spin-spin correlations obtained with QMC to calculate the staggered /-orbital magnetiza- 
tion in the thermodynamic limit. When this quantity is non-zero, the ground state of the 
PAM has long-range antiferromagnetic order. In addition, we have evaluated non-equal time 
correlations to extract the dynamic /-orbital spin and charge susceptibilities as well as the 
low temperature dependence of the total uniform spin and charge susceptibilities in order 
to determine the spin and charge gaps. For an antiferromagnetic Mott insulator the charge 
gap is finite but there are gapless Goldstone antiferromagnetic spin excitations. However, 
in a Kondo insulator there are gaps for all excitations (with the charge gap larger than the 
spin gap) and consequently no long-range order. In the PAM we expect the small exchange 
interaction phase to be well described by a Mott insulator with antiferromagnetic order, 
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while the larger exchange interaction phase should be well described by a Kondo insulator 
with spin-liquid behavior. In addition, in the limit of a large exchange interaction there is 
a cross-over to a band-insulator regime where the spin and charge gaps become identical. 

At finite temperature we find two different behaviors depending on whether the parame- 
ters are such that the ground state is in the order or disordered part of the phase diagram. For 
small values of the exchange coupling we find a metallic behavior, with a Kondo resonance 
developing in the /-band density of states for temperatures below the Kondo temperature 
Tk- This resonance is then suppressed at a lower temperature, leading to insulating behav- 
ior, by the formation of a gap when the antiferromagnetic correlations reach the size of the 
system. The staggered spin susceptibility tends to diverge at the same temperature at which 
the charge gap forms in N^{uj). The uniform spin susceptibility displays a maximum below 
T K and tends to saturate at a finite value as T goes to zero. On the other hand, the uniform 
charge susceptibility also displays a maximum but then decays quickly to zero, indicating 
the existence of a charge gap. For large values of the exchange coupling the /-band density 
of states shows an insulating gap for all finite temperatures and the absence of antiferro- 
magnetic order. Both the uniform q — > spin and charge susceptibilities decay to zero as 
the temperature is lowered, indicating the presence of a gap. 

In order to determine the Uf-V boundary separating the antiferromagnetically ordered 
from the disordered phase we have calculated the /-band spin correlation function 



as functions of temperature and system size N. In particular, for a given lattice size, we 
found that Sf(ir, ir) saturates for large values of (3. We have used the saturated values along 
with the spin-wave scaling forms 




(2) 



with ml = nL — nf,, and the /-band antiferromagnetic structure factor 



s f (7r,n)=j:(-iy* +iy c f (i x j y ), 



(3) 




(4) 
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and 

(5) 



mj = 3c f (y/N/2, y/N/2) + O (^=j > 



to extract the staggered magnetization mj of the /-band. Figure 1 shows the results of 
these calculations for two values of the hybridization V = 1.0 and V = 1.2 with Uf — 4. 
The system sizes ranged from 4 x 4 to 8 x 8 lattices. For V = 1.0, the intercept of a linear 
fit to the points gives a finite value for mj in the infinite system. However, for V = 1.2, 
S-^ (71,71) does not appear to scale with system size N and, as shown in Fig. 1, we conclude 
that m'j = 0. Similar conclusions can be reached by scaling c^(\fN /2, /2) according to 
the spin- wave theory result. In Fig. 1 we show that the extrapolated value for mj obtained 
from Eq.(||) is consistent with the value obtained from S^(tt,7t) for V = 1.0. Again, for 
V = 1.2 it appears that mj = 0. 

We have carried out a similar scaling analysis for several values of Uf and in Fig. 2 we 
show the Uf — V phase diagram obtained in this manner. The boundary is shown as a solid 
line passing through the three points obtained from the QMC data at Uf = 2, 4, and 6. 
For small values of J ~ V 2 /Uf we have carried out a mean- field calculation assuming that 
{ n ia) = \ + cr^(—iy. The dotted curve in Fig.2 shows where the mean-field mf vanishes. 
The antiferromagnetic region obtained from MFT is significantly larger than that found 
from QMC due to the effect of fluctuations which are neglected in MFT. 

In order to extract dynamic properties of the system, we have also calculated imaginary- 
time properties. Using a maximum entropy procedure [|14]] to analytically continue the QMC 
data, we have determined the density of states and the imaginary part of the spin and charge 
dynamic susceptibilities. We have inverted the integral relation relating the single particle 
Green's function G{j(r) = — (T t q(t)c](0)) to the density of states N^(u), 
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In Fig. 3 we show N^(u) for several values of (5 for Uf = 4 and V = 0.75. As the temperature 
is lowered from (3 = 4 to f3 = 12 we see a Kondo resonance peak developing at the Fermi level, 



corresponding to the screening of the local moments by the conduction electrons. However, 
for finite values of V, as the temperature is further lowered, we find that the Kondo peak is 
suppressed until at f3 = 20 a well developed gap is formed. On the same lattice size and for 
the same range of temperatures we find that the antiferromagnetic structure factor grows 
rapidly before saturating to a finite value at around f3 — 20. Thus, we associate the gap 
with the long-range antiferromagnetic order which sets in at the same temperature for this 
lattice size. For larger values of V the behavior of N^(uo) is quite different, since we do 
not find any resonance peak at uj = but rather an insulating gap is always present at all 
temperatures. It is found that this gap grows with V for any fixed Uf. 

We have also analytically continued the imaginary-time /-orbital spin and charge corre- 
lation functions XiA T ) = (T T Tn{ (T)m{ (0)) and Ilf^r) = (T T n{(r)n{ (0)) in order to calculate 
the imaginary part of the dynamic spin and charge susceptibilities Im%^(a;) and lmU^(uj), 
respectively. It can be shown that Imx^co) is related to Xu{ T ) by 

with a similar expression relating ImU^cu) to Hf^r). In Fig. 4(a) we show -Im%^(a;) for 
several choices of V on a 6 x 6 lattice. The spin gap only opens for values of V above the 
critical point, where there is no antiferromagnetic long-range order. In Fig. 4(b) we show 
-lmW(uj) with the same parameters as in Fig. 4(a). In this case it is clear that the charge 
gap A c is present also for values of V below the critical point. In the inset we report the 
value of the ratio of the spin gap to the charge gap A s /A c as a function of V. We find that 
A c > A s for all values of V. However, for the larger values of V we can clearly see that 
the two gaps tend to be equal. This behavior is an indication that the system crosses into 
a band-insulator regime. 

Evidence for charge and spin gaps can also be seen in the temperature dependence of 
the uniform charge and spin susceptibilities. The spin susceptibility is defined by 

Xtot(? = 0)= [ P dr(m(r)m((})}, (8) 
Jo 
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and the charge susceptibility is similarly defined by 



ntot(? = 0) = / dr [(n(r)n(O)) - (n) 2 } . (9) 
Jo 1 J 



Here m = ^ J2i (nq - and n = Ei (^t + witn n ^ = n »C + n t- In Fig. 5 we 
show xtot an( i n^ot versus temperature for several values of the hybridization energy V at 
a fixed Uf — 4. For smaller values of V we find that n^ t becomes very small at a finite 
temperature, while xtot peaks at a temperature close to the one at which we find the Kondo 
resonance in N^(uj) and then approaches a finite value at T = 0. This behavior is consistent 
with the fact that for V — 1.0, ImU^uj) has a gap but Imx^(cj) is gapless. 

In conclusion, we have shown that for small values of V 2 /Ut the ground state of the PAM 
is an insulator with long-range antiferromagnetic order characterized by a finite charge gap 
and gapless spin excitations. As V 2 /Ut increases the long-range order is destroyed and the 
system exhibits a spin-liquid behavior. The spin-liquid state is characterized by both a spin 
gap and a charge gap with A c > A s . When the hybridization V increases further and V 2 /Ut 
becomes large, the system crosses over to a band-insulating state in which A s approaches 
A c . Work is currently in progress to determine the behavior of the doped system. 
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FIGURES 



FIG. 1. The square of the staggered magnetization m'j for various y/N x yN lattices as 
a function of 1/ y/~N. The solid symbols are extracted using Eq. (^) and the open symbols 
using Eq.(|5p. The circles are for V = 1.0 and the triangles for V = 1.2. Here Uf — 4 and 
(3 = 20. For comparison, the solid square indicates the value of mj for the 2D Hubbard 
model with U = 4. 

FIG. 2. The U f-V phase diagram of the half-filled PAM at T = 0, showing the boundary 
between the antiferromagnetic (AF) and spin-liquid (SL) phases. The solid line is a guide 
to the eye through the QMC points. The dotted line is the MFT boundary. 

FIG. 3. The /-band density of states N*((jj) on a 6 x 6 lattice with Uf — 4 and V = 0.75 
for several temperatures. The dashed line is for j3 — 4, the dashed-dotted line for /3 — 12 
and the solid line for j3 = 20. 

FIG. 4. The imaginary part of (a) the spin susceptibility -Imx (u) and (b) the charge 
susceptibility -Indira;) for f3 = 16 on a 6 x 6 lattice with Uf = 4. The ratio of the spin and 
charge gaps A s /A c versus V is shown in the inset. 

FIG. 5. The uniform (a) spin susceptibility Xtot(<? = 0) an d (b) charge susceptibility 
ntot(<? = 0) as a function of temperature T on a 6 x 6 lattice with Uf = 4. 



10 



o 






o 


LO 


o 


LO 


o 


LO 


o 


CO 


CO 


CM 


C\2 






o 


o 


• 

o 


• 

o 


• 

o 


• 

o 


• 

o 


• 

o 


• 

o 


• 

o 



( ro )jN 



